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Abstract 



We have obtained a set of coupled differential equations from the continuous limit of the transfer 
matrix method. Decoupling such a set of equations yields an extension to the Wentzel-Kramers- 



Brillouin (WKB) approximation for the Schrodinger equation with the position-dependent effective 



mass (PDEM). In the classically allowed region, the decoupling is to ignore the reflection resulting 
from the variations of both the potential and effective mass. By considering an infinite- well example 
with the PDEM, it is shown that the extended WKB approximation can provide not only an 
estimation to eigenenergies, but also an analytic form to approximate wavefunctions. 
* Corresponding authors: sdchao@spring.iam.ntu.edu.tw (S. D. Chao); drhang@mail.nsysu.edu.tw 
(D. R. Hang) 
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I. INTRODUCTION 



The one-dimensional Schrodinger equation with the position-dependent effective mass 
(PDEM) has been introduced for many microstructures of current interest. To 
understand such an equation analytically, exactly solvable examples have been discussed in 



the literature. 



, |6|, |7|] On the other hand, the transfer matrix method {si, ^, [h], ll|, 12 1 is a 
powerful numerical way to calculate the eigenvalues and eigenfunctions of the Schrodinger 
equation with the PDEM. To perform this numerical method, the space is divided into a 
series of slabs so that the effective mass and potential can be approximated as constants in 
a slab. By considering the conservation of the flux, it has been suggested that we shall use 
the BenDaniel condition to match an eigenfunction at each boundary between two adjacent 
slabs. 0, [9I, [h], Q In Ref. 1, Ou, Cao, and Shen derived a dispersion relation for the 
one-dimensional Schrodinger equation with the PDEM by considering the analytic transfer 
matrix method. Such a dispersion relation contains an integral term of Wentzel-Kramers- 
Brillouin (WKB) type. It is well-known that WKB approximation 



17l | provides 



important classical interpretations in the developement of quantum mechanics. Hence it is 
interesting to probe the validity of WKB approxmation after incorporating the variation of 
the effective mass. 

In this paper, a set of coupled differential equation is obtained by considering the continu- 
ous limit of the transfer matrix method. Decoupling such a set of equations, the approximate 
wavefunction of the analytic form 



ip(x) 



'm*(x) 
k(x) 



exp[±iS(x)/h] 



(1) 



with 



S(x) = J k{x')dx' (2) 

is derived in section II when the effective mass m*(x) depends on the position x. Here V(x) 
is the potential, Ti equals Planck constant divided by 27r, k(x) = y^2m*(x)(E — V(x))/h, 
and E is the eigenenergy. The phase S(x) is of the same form as that in the well-known 
WKB approximation except that m*(x) is position-dependent. Therefore, we can extend 
such an approximation for the one-dimensional Schrodinger equation with the PDEM from 
the transfer matrix method. When the effective mass becomes constant in space, the ap- 
proximate solution given by Eq. (1) is reduced to the conventional WKB-type function. 



As discussed in section III, the decoupling in the classically allowed region is to ignore the 
reflection due to the variations on the effective mass and/or potential. From our study, 
it becomes clear how to relate the numerical solution obtained from the transfer matrix 
method to the analytic WKB approximation. By considering an infinite-well example, it 
is shown that the extended WKB approximation can provide not only an estimation to 
eigenenergies, but also an analytic way to understand the eigenfunctions when the effective 
mass is position-dependent. The conclusion is made in section IV. 



II. TRANSFER MATRIX METHOD AND APPROXIMATE WAVEFUNCTION 
OF WKB-TYPE 



When the effective mass m* depends on the position x, it has been shown that the 
one-dimensional Schrodinger equation should be modified as 



d hf 



d 



ip(x) + V(x)if)(x) = Eip(x). 



(3) 



dx 2m* (x) dx 

In this paper, we assume that m*(x) > 0. In addition, let m*(x) and V(x) be analytic 
functions. For convenience, in this section we only consider the classically allowed region, 
where 



E - V(x) > 

at any point x. To extend WKB approximation, the equation 



(4) 



d 


t(x) 




t(x) 




= T(x) 




dx 


r(x) 


r(x) 



(5) 



will be derived in this section by considering the continuous limit of the transfer matrix 
method. Here the position-dependent matrix 



T(x) = 



-ixk'(x 

m*(x) I k(x) 



m* (x) ( k(x) \ m*(x) ( k(x) ^ p -2ixk(x) 

2k(x) \m*(x)j ' 2k(x) \m*(x) 
I 



D 2ixk(x) 



ixk'(x) — 



m*(x) I k(x) 



( , (6) 

2k(x) y m *(x)J ° ''' *^ru y^j ^ m *( x ) 

t(x) and r(x) are two complex functions, and in this paper the notation f'(x) presents the 
derivative of any function f(x). With some calculations, we can see that Eq. (5) is equivalent 
to Eq. (3) if we set 



ij){x) = t(x)e ik{x)x + r(x)e~ ik{x)x . 



(7) 



Then approximate solutions of WKB type are obtained by decoupling Eq. (5). Therefore, 
we can relate WKB-type solution to the transfer matrix method. 

The transfer matrix method is a powerful numerical approach to solve Eq. (3). 
To perform such a method, as shown in Fig. 1, we can divide the space into a series of 
thin slabs so that both the effective mass and potenital can be approximated as constants 
in each slab. For convenience, each slab is labelled by an integer j sequently from the left to 
right. In the classical allowed region, the wavefunction ip(x) in a slab can be approximated 
as the linear superposition of two plane waves. Let xj be the center point of j-th slab, and 
approximate the effective mass, potential, and wave number as m* = m*(xj), Vj = V(xj), 
and kj = k(xj) in such a lab, respectively. Then the wave function is taken as 



ip(x) 



tje ikjX + rje ikjX 



(8) 



in the j-th slab. The coefficients tj and rj are the traveling components for two different 



directions. We shall consider the BenDaniel condition [13] to relate the coefficients of adja- 
cent slabs. Q] Based on such a condition, as shown in Appendix, the coefficients tj+i 
and Tj + i are related to tj and rj by 



with the transfer matrix 



75 



tj+i 



l/i _ ^ m *j+i \ i{kj+k j+1 ) yj 1(1 , kj -j(fc,— k j+ i)y 

m* k j+1 ) C ' 2^ ^ m* k J+1 > C 



1(1 _ ^ m J+1 
2 V m* kj +1 

k 3 _^+l- 



-i(kj+kj +1 )yj 



(9) 



(10) 



after some calculations. Here t/j is the point separating the j-th and (j + l)-th slabs, so 
Xj+i = {yj + yj+i)/2. It is shown in Appendix that we can relate the transfer matrix Tj to 
T(x) by 



% = 1 + ^)^+0(^1 



(11) 



with Axj = Xj + \ — Xj and 



1, 



0, 1 



(12) 
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Let Atj =tj + \ — tj and Arj 
of equation, 



'3+ 



i — Tj. From Eqs. (9) and (11), we have the following set 



At, 








= nvs) 


Xj 


_A rj _ 


r x- 

Xj 



+ o(A Xj ). 



(13) 















t 


>i) 




<—> 






_ T 3 _ 




r 


[ x j) _ 



By shriking the widths of slabs so that Axj approaches zero, the above equation can be 
reduced to Eq. (5) with the following correspondence 

(14) 

Thus we can obtain Eq. (5) by considering the continuous limit of transfer matrix method. 
From Eqs. (8) and (14), we can see why the wave function is determined by Eq. (7) after 
obtaining r(x) and t(x). 

In Eq. (5), the coefficients t(x) and r(x) are coupled to each other by the off-diagonal 
terms of T(x). With some calculations, we can see from Eq. (11) that these terms comes 
from the off-diagoal terms of Tj, which corresponds to the reflection due to the variations 
on potential and/or effective mass , [loj]. If we ignore the reflection by neglecting the 
coupling bwtween t(x) and r(x), we have the following decoupled equations 



—t(x) ~ \-ixk'(x) - ^ ( 

dx [ ' [ [ ' 2k{x) \m*(x) 

dx [ ' 1 1 ' 2k(x) \m*(x) 



]t(x) 
r(x). 



Solving the above two first-order differential equations, we have 

t f x )jKx)x 



' m*(x) 



k(x) 



exp(iS(x)/h) 



r{.v)(i ~ lk{ - x > ~ c 2 J^^exp(-iS(x)/h). 



k(x) 



(15) 
(16) 

(17) 
(18) 



Here c\ and c 2 are two constants. From Eq. (7), the approximate solution is of the form 

(19) 



' TTl* [ X) TYl* ( X) 

ijj(x) ~ c\ a I — ——— exp(iS(x)/h) + C2\ exp(—iS(x)/h). 



k(x) y k(x) 

We shall set ci=0 (c2=0) such that t(x)=0 (r(x)=0) for the traveling wave moving to right 
(left), and obtain the WKB-type function given by Eq. (1). Therefore, the analytic func- 
tions of WKB-type can be related to the numerical transfer matrix method by ignoring the 
reflection in the continuous limit. 



5 



III. DISCUSSION 



In the last section, the approximate solution of WKB-type is derived for the classically 
allowed region by ignoring the off-diagonal terms of T(x) to decouple Eq. (5). With in- 
creasing eigenenergy E, the ignorance is reasonable as k(x) » 1 because the off-diagonal 
terms contain the factor e ±2lk ( x ) x ; which oscillates quickly under large k(x). On the other 
hand, the diagonal terms of T(x) do not contain such an oscillating factor. Hence Eq. (1) 
provides a good approximation for the states with high energies as in the conventional WKB 
approximation. Because the transfer matrix Tj can be related to T(x) by Eq. (11), we can 
expect that the numerical solution obtained by the transfer matrix method can be reduced 
to the WKB-type function as k(x) » 1. 

It has been shown that the dispersion relation of WKB-type can be used to estimate the 
eigenenergies of Schrodinger equation with the PDEM. |l|] To further probe the extension 
of the WKB approximation to systems with PDEM, we consider the infinite quantum well 
where 

m 1 -m 2 mi+m 2 

m (x) = x J \ 20 

K J 2a 2 v ' 

and V(x) — if \x\ < a while V(x) = oo if \x\ > a. We do not need to consider how the 

phase changes at the turning points in this case, 

HQ 

and just need to set eigenf unctions 
as zeros at x = ±a. Hence we shall take c\ = —c 2 = c in Eq. (19) to obtain the approximate 
eigenfunctions of WKB type 



Mx) ~ cJ^f^-sin(iS(x)/n) (21) 

V k ^ x > 

if we set S(x) = J^ a \j2m*(x')(E — V(x'))dx' and require S(a) = nrc. Here n is a nonnega- 
tive integer. Then we can obtain the WKB dispersion relation 

E n = ^ with the integer n = 1, 2, 3, ... 22 

128(mi /2 - mf) 2 a 2 

On the other hand, we can reduce Eq. (3) as 

dy 

in such an example if we set 

(m 1 -m 2 )E 1/3 m 1 + m 2 . 

V — —( 9 ) ' (x H a) 

V h 2 a ' y mi -ma ' 



2 u(y) = yu{y) (23) 



and u(y) = ij)'(x)/m*(x). The solution of the above equation is the linear combination of Airy 



functions, 



17 



20J and thus it is easy to obtain the exact eigenvalues and wavef unctions. In 



Table 1, we compare the exact and approximate eigenenergies when m x = 0.1 m , m 2 = 0.2 
m , and a = 5 nm. Here uiq is the rest mass of electrons. We can see from table 1 that the 
errors of WKB approximation are reduced with increasing n. The error for the n=4 excited 
eigenlevel is only about 0.14%, and we can see from Fig. 2 (a) that the corresponding WKB 
wavefunction is very close to the exact one. (In Fig. 2, both the approximate and exact 
wavefunctions are normalized such that f_ \ip(x)\ 2 dx = 1.) Therefore, the extended WKB 
approximation provides not only a good estimation to the eigenenergies, but also a good 
way to understand the wavefunctions. Figure 2 (b) shows the square of n = 4 wavefunction. 



The oscillating amplitudes, in fact, are proportional to the factor y / m*(x)/k(x) as expected 
from Eq. (21). 

For the tunneling through a classically forbidden region where E — V(x) < 0, it is known 
that we can obtain an estimation for the tunneling probability from WKB approximation 
if the (effective) mass is constant. With some calculations, we can obtain Eq. (5) with 
the same r(x) for the classically forbidden region from the continuous limit of the transfer 
matrix method. (It should be noted that k(x) becomes imaginary.) In addition, the WKB- 
type approximate solution given by Eq. (1) can still be obtained by decoupling Eq. (5). In 
fact, it is shown in Ref. 6] that Eq. (3) can be transformed to the conventional Schrodinger 
equation 

l£- 2 -^(v(y)-EMy) = v (24) 



if we set y = J ^m*{x)dx ) (f)(y) = i/j(x)/m* 1 ^' l (x), and V(y) = V(x) + F(x). Here 

/ x h 2 r^-2 m *(x) 7 / -fm*(x),„ 1 
p(x) = \ dx ( dx ) } 

8m*(x) m*(x) 4 m*(x) 

With some calculations, in fact, we can see that Eq. (1) can also be obtained from the 
conventional WKB-type function by considering such a transform after ignoring F(x). The 
ignorance is reasonable when m*(x) varies so slowly that its derivatives are small. Because 
WKB approximation may provide a suitable approximation to the above equation in the 
classically forbidden region, Eq. (1) also provides an extension to the WKB approximation 
when E - V(x) < 0. 
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Just as in the conventional WKB approximation, the reflection is neglected if we use 
Eq. (1) to approximate the wavef unction. If the effective mass or potential has sharp 
jumps at some interfaces, we cannot ignore the reflection and Eq. (1) does not give us a 
good approximate solution in the whole space. But we may apply Eq. (19) piecewisely 
to approximate the wavefunction. For example, consider a finite-well problem where the 
potential and effective mass have sharp jumps at two interfaces A\ and A 2 , as shown in Fig. 
3. We can divide the space into three regions denoted by (I), (II), and (III), respectively, and 
approximate the wavefunction by Eq. (19) in each region. Then we just need to consider the 
BenDaniel condition at A\ and A 2 to obtain the approximate solution in the whole space. 

IV. CONCLUSION 

In this paper, we derive a set of first-order differential equations corresponding to the 
continuous limit of the transfer matrix method. Decoupling such a set of equations, an 
approximate wavefunction of WKB-type is obtained for the one- dimensional Schrodinger 
equation with the position-dependent effective mass. In a classically allowed region, the 
decoupling is to ignore the reflection induced by the variations of the effective mass and 
potential. Such ignorance is reasonable in the high-energy limit. From our derivation, it is 
clear how to relate WKB approximation to the numerical solution obtained by the transfer 
matrix method when the effective mass is position-dependent. By considering an infinite- well 
example, it is shown that such an approximation can provide not only a simple estimation 
to eigenenergies, but also an analytic form to approximate the eigenfunctions. 
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Appendix 



In the transfer matrix method, the BenDaniel condition is taken into account to match 
ip and ip'/m*, the wavefunction and its derivative divided by the effective mass, at the 
boundary of each slab. At the boundary yj separating the j-th and (j + I)-th slabs, we have 

t j+1 e ik * +iy i + r j+1 e- ik i +iy i = t j e ik M + r je ik ^ (25) 

by matching ip and 

ih+L t e ik j+m _ ih+L r e -ik j+m = % ^L t ^m _ i^L T .^m ( 2 6) 

m* j+1 ]+ m* j+1 J+L m* J m* 3 y ' 

by matching ip'/m* if ip is approximated by Eq. (8). We can obtain Eqs. (9) and (10) by 
relating tj +1 and r J+ i to tj and rj from the above two equations. 
By introducing the function h(x) = k(x)/m*(x), the factor 

ki m ^ = = 1 - ^Ax + o(Ax 2 ). (27) 



m* k j+1 h(x j+1 ) h(yj) 

Based on the above equation, we can obtain Eq. (11) as the first-order approximation with 
respect to Ax since e ±l( - kj+1 ~ kj ^ yj ~ 1 ±iyjk'(yj)Ax and e ±l( - kj+1+kj ^ Vj ~ e ±2tk( - yj ^ yj . 
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Figure caption 

Fig.l To perform the transfer matrix method, the space is divided into a series of slabs, in 
each of which both the potential and effective mass are approximated as constants. The 
points Xj and Xj+i present the centers of the j-th and (j + l)-th slabs, and yj is the point to 
separate these two slabs. 

Fig. 2 (a) The red dash line and black dot line corresponds to the exact wavefunction and 
WKB approximate one. (b) The solid line is the square of the wavefunction and the dot 
line is the envelope function proportional to m*(x) 1 ^ 2 . 

Fig. 3 A finite-well problem where both the effective mass and potential have sharp jumps 
at the interfaces Ai and A 2 . The regions denoted by (I), (II) and (III) are divided by these 
two interfaces. 

Table caption 

Table 1 The comparison of WKB and exact eigenvalues for the solvable example discussed 
in the text. 
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n 


WKB(eV) 


Exact (eV) 


Error (%) 


1 


0.0253 


0.0258 


1.93 


2 


0.1012 


0.1018 


0.55 


3 


0.2278 


0.2283 


0.25 


4 


0.4049 


0.4055 


0.14 


5 


0.6327 


0.6333 


0.09 


6 


0.9111 


0.9117 


0.06 


7 


1.2401 


1.2407 


0.05 


8 


1.6197 


1.6203 


0.04 


9 


2.0499 


2.0505 


0.03 


10 


2.5308 


2.5313 


0.02 



Table 1 C. F. Huang et al. 
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Fig. 2 (a) C. F. Huang et al. 
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Fig. 3 C. F. Huang et al. 



